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ABSTRACT 

The luminosity function of active galactic nuclei has been measured down to luminosities ~ 10 42 ergs s" 1 
in the soft and hard X-rays. Some fraction of this activity is associated with the accretion of the material 
liberated by the tidal disruption of stars by massive black holes. We estimate the contribution to the X-ray 
luminosity function from the tidal disruption process. While the contribution depends on a number poorly 
known parameters, it appears that it can account for the majority of X-ray selected AGN with soft or hard X- 
ray luminosities < 10 — 10 44 ergs s" 1 . If this is correct, a picture emerges in which a significant portion of the 
X-ray luminosity function of AGN is comprised of sources powered by tidal-disruption at the faint end, while 
the sources at the bright end are powered by non-stellar accretion. Black holes with masses <2x 10 6 M Q could 
have acquired most of their present mass by an accretion of tidal debris. In view of the considerable theoretical 
uncertainty concerning the detailed shape of the light curves of tidal disruption events, we focus on power-law 
luminosity decay (as identified in candidate tidal disruption events), but we also discuss constant accretion rate 
models. 

Subject headings: accretion — black holes — cosmology: observations — galaxies: active — galaxies: nuclei 
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1. INTRODUCTION 

A main-sequence or giant star passing near a massive black 
hole ( MBH) of mass M hh < 10 8 M© is disrupted by the tidal 
forces (iHillsl 19751) . The disruption and the subsequent accre- 
tion of stellar debris by the black h ole are potential power- 
ful sources of X-ra y emi ssion ( e.g.. | Meszaros & Silklll977t 
Young et alJ Il977t IReesI 119881 I199oTk The emission can 
be produced in various phases of tidal disruption, including 
the initial compression of the star ( Kobavashi et al. 2004), 
the interaction of ejected debris with the ambient medium 
JKhokhlov & Merial l'l996). in tidal stream collisions and de- 
bris fallback onto a nascent accretion disk jKochanelj l 19941 
iLee et alJH996t iKim et aljfl999t iLi et al.ll2002l). and in the 
accre t ion of disk m aterial onto the black hole llCannizzo et ail 
1990: IUlmerl [l999). The purpose of the present study is to 
estimate the contribution of activity associated with tidal dis- 
ruptions to the luminosities of typical active galaxies. 

Recently, surveys with HEAO 1, ASCA, Chandra, and 
XMM-Newton have yielded luminosity functions (LFs) of 
active galaxies at low and high redshi ft in the hard 
(e.g.. lUeda et al] 120031 iBar ger et all 12005) and soft (e.g., 
Hasinger et al. 2005) X-rays. The measured X-ray LFs extend 
to lower luminosities than their optical counterparts because 
contamination with host galaxy light makes optical detection 
of low-luminosity AGN difficult. The LFs are approximately 
broken power laws dN / d In Lx ocL^(Lx+ix*) ly_A ' with vfa—2 
and -1 < fj, < and Lx* ~ 10 43 - 10 44 ergs s" 1 . The power- 
law behavior at the bright end has been attributed to the in- 
trinsic structure of quasar luminosity histories. The luminos- 
ity histories (light curves) are convolve d with the black hole 
mass function to calc ulate the LF (e.g., Wvi the & Loebl2003l: 
Hopk ins et al.l l2005). The origin of the LF at the faint end, 
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however, remains unknown. 

At least a fraction of low-luminosity AGN must be powered 
by stellar tidal disruption. Here we estimate the contribution 
of these sources to the X-ray LF of AGN. This estimate is an 
essential step to understanding the origin of activity in AGN 
with black hole masses < 10 7 M© and ultimately to under- 
standing the origin of low-mass MBHs. 

Theoretical investigations of the stellar tidal disruption 
process, and the detection X-r ay flares that are candi- 
date tidal disruption e vents (e.g., iKomossa & GreinerlH999t 
iKomossa et al.ll2004bt lHalpern et alJl20 04). suggest that the 
nuclei in which stellar tidal disruption has occurred remain 
luminous for 1-10 yrs after disruption or longer. The inte- 
grated luminosity of the AGN depends on the mass of the 
stellar debris that remains bound to the black hole, and on 
the fraction of this mass that ends up accreting onto the black 
hole in a radiatively-efficient fashion. The bound mass can be 
a frac tion f ~ 0.25-5 of the initial stellar mass jAval et alJ 
2000). 

The bound mass inferred in the giant X-ray flare in NGC 
5905 using a bl ack-body fallb ack model is much smaller than 
the solar mass tLi et al.l l2002): in two other candidate flares, 
RX J1624.9+7554 and RX J1242. 6-1119A, the light c urves 
are consistent with masses ~ Mq. iDonlev et alJ (120021) con- 
strained the rate of outbursts by comparing the ROSAT All- 
Sky Survey with archival pointed observations. They inferred 
a rate of ~ 10~ 5 yr" 1 per galaxy. We adopt a theoretical es- 
timate for the tidal disruption rate that is significantly higher 
(eq.D- 

T he ou tline of our calculation of the AGN LF is as follows. 
In § 12.11 we estimate th e ma ss function of MBHs and the rate 
of tidal disruptions. In § 12.21 we present sim ple m odels for the 
light curves of tidally disrupted stars. In § 12.31 we calculate 
the LF and explore its dependence on various parameters in 
the model. In §[3] we discuss implications of our estimate for 
the understanding of low-luminosity, X-ray selected AGN. 

2. CALCULATION OF THE LUMINOSITY FUNCTION 
2.1. Tidal Disruption Rate 
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Stars passing within distance (r/ z Mtf l /m ir ) 1 ' 3 R- k of the black 
hole, where is the mass of the star, R+ is its radius, and 
rj ~ 1 is a numerical factor, are tidally disrupted. Here we 
estimate the tidal disruption rate in a galaxy with black hole 
of mass Mbh and stellar velocity dispersion a. The disruption 
rate depends on the detailed structure of the stellar nucleus 
of the galaxy. The galaxies relevant to this study are those 
with M bh < M max = lO 8 M Q (/^/fl Q ) 3/2 (m,/M )- 1/2 which 
can disrupt stars. 

In spheroids as faint as that of the Milky Way, the relax- 
ation time within the black hole's radius of dynamical in- 
fluence rbh ~ GM bb /er 2 is short enough that a collisionally- 
relaxed distribution, p oc r~ 7 , where 7 ~ 1.5— 1 .75, i s 
set up jBahcall& Wol3ll976l Il977t iMerritt & Szel 120051) . 
This is consistent with what is seen a t the Galactic Center 
JGenzel et alJ2003t l Schodel et alJ20 06). But just beyond r bh , 
the Galactic Center density profile steepens to p oc r~ 2 , and 
this is also the slope observed with the Hubble Space Tele- 
sco pe near the centers of all but the most l uminous galax- 
ies (Gebha rdt et alJl996tlFerrarese et alJ2006l) . Furthermore, 
nuclear relaxation times in gala xies with M bb > W Mq gen- 
erally exceed a Hubble time (Faberetal. 1997) and so a 
Bahcall-Wolf cusp would not form. 

The rate of stellar disruptions in stellar den sity cusps has 
been s tudied at various levels of detail ("e.g.. iFrank & Reesl 
[l976t iLiehtman & Shapiro! 119771 iCohn & Kulsrudl 119781 
Magorrian & Tremaine 1999). We adop t the most recent esti- 
mate of the rate in a p oc r~ 2 cusp J Wang & Merritl2004l) 
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where m + and R+ are the mass and radius of the tidally dis- 
rupted stars. 

The disruption rate in equation Q must be weighted by 
the stellar mass function in the spheroid at distances from 
the black hole corresponding to the initial pericenter radii of 
tidally disrupte d stars. These r adii are typically similar to r b h 
(Magorrian & Tremaine 1999). If the stellar m ass function re - 
sembles the "universal" initial mass function (Kroupa 2001), 
the disruption is dominated by ~ O.1M stars. However, 
the mass function near the black hole may differ from the uni- 
versal initial mass function because higher-mass stars dynam- 
ically segregate closer to MBH, and because stars with an un- 
usual initi al mass function may f orm in situ near the MBH 
(e.g.. lMilosavlievic & Loebl20 04. and references therein). 

Mounting observational evidence indicates that the black 
hole mass and the velocity dispersion of the host stellar 
spheroid are tightly correlated with best-fit relation M b h « 
1.7 x 10 8 M Q (a/200 km s" 1 ) 4 - 9 fe.g.. lFerrarese & Fordll20l)5l 
and references therein). We explore departures from this fidu- 
cial relation by parameterizing the logarithmic slope of the 
relation M b h oc a" . 

No direct constraints on the cosmic density of MBH with 
masses < 10 7 Mq exist. The majority of these MBH should 
be in spiral galaxies. The density can be estimated by as- 
suming a relation between the bulge luminosity and the black 
hole m ass of the form L bu i ge = AM^ h . Following Fe rraresd 
(2002), we substitute this in the Schechter LF <&{L)dL = 
^{LlU) a e~ Llu dLjU to obtain a cosmic mass function of 
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black holes 

y{M hh )dM hh = * (^) 

where f = k$ , M+ = (1 .27 0U/A) l l k , and = L hn w./L ~ 
0.3 for the Hubble type Sab ( Simi en & de Vaucouleursl 19861) 
which we use as reference. The latter is justified be- 
cause spiral gal axies dominate galaxy LF below ~ (e.g., 
Nakamura et al 1l200l . For the relation between the blue lu- 
minosity of the galaxy a nd the black hole mass we adopt 
L = 3900L Q (M bh /M Q )° 79 ( iMarconi & Huni2003l) . 

The number density of ~ 10 8 M Q black holes esti- 
mated using equatio n agrees with a similar estimate in 
Marc oni et alJ ll2004l) . while the number density of ~ 1O 6 M 
black holes exceeds their estimate by a factor — 3-4. This 
may be a consequenc e of our using the ff -band galaxy LF 
as reference, whereas Marcon i et alJ ( 120041) use the A'-band 
LF that is flatter at the faint end. However our adoption of a 
steeper M b h-<r relation than they do has an opposite effect. 

2.2. The Light Curve of a Tidal Disruption Event 

The light curves o f candidate tidal disruption events 
( iKomossa et alJl2004bl) are characterized by rapid fall in X- 
ray luminosity on time scales of months to a year, followed by 
a gradual continued fading over the following decade. Three 
of the four candidate events (NGC 5905, RX J 1420+53, RX 
J 1 242— 11) are consistent with luminosity decay Lfb oc t~ 5 / 3 
expected if the X-ray emission is produced during the fallback 
of stellar debris onto a na scent accretion disk ( Phinne\ fl989t 
lEvans & Kochanekll989l) . The fourth event RX J 1624+75 ex- 
hibits what appears to be faster initial decay followed by very 
slow fading. In all cases, the total drop in X-ray luminosity is 
by a factor - 10 2 -10 3 over about a decade. 

The above candidate tidal disruption events have been se- 
lected for characteristic, "outburst"-like light curves (fast rise 
and slow decay). However the true luminosity evolution of 
tidal disruption events is unknown. If the emission during the 
fallback stage is absorbed or otherwise suppressed, the X-ray 
light curve will be dominated by thin disk accretion and a 
quasi-steady flux over a longer period will be expected. Such 
events would remain undetected in existing searches that tar- 
get outburst activity. To maintain generality, therefore, we 
consider fallback-like power-law light curves, and also dis- 
cuss constant accretion rate models. 

The power irradiated during fallback was derived in lLi et all 
J2002I) 
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where the tidal disruption occurs at t = 0, and the luminos- 
ity reaches peak value and starts decaying as in equation (0 
at time f pea k ~ O.O2(M bh /lO 6 M Q ) 1 /2( m ,/M )- 1 (^/« ) 3 / 2 yr 
after tidal disruption. 

Although equation was derived for the fallback process, 
we interpret it as a representative power-law decay model for 
the luminosity evolution, independent of the emission mech- 
anism To test the sensitivity to the power law index of lumi- 
nosity decay, we consider decay of the form L bo i oc f~ A , where 
A > 1 is a constant. We also test the sensitivity to the value of 
fpeak by scaling it via f pea k(0 = £ f peak(0). In both cases we keep 
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the total energy irradiated between t = f pea k and ? = oo constant 
and equal to /ra*c 2 . 

We also discuss an alternative to power-law light curves. 
Accurate theoretical modeling of a time-dependent thin disk 
emission is difficult. We therefore adopt a crude quasi-steady 
accretion scenario, in which we assume that accretion at con- 
stant rate proceeds "while the supplies last," i.e., until all the 
bound debris has been accreted. This is inevitably an over- 
simplification. Taking into accoun t the redistribution of mass 
in the disk, Cannizzo et al. ( 1990) estimate that under certain 
assumptions the light curve will be Ldisk °c t~ 1 - 2 . On the other 
hand, if outer parts of the disk become thermally un stable 
and n eutral, the accretion may be delayed ( Menou & Ouataertl 
12001 . 

In the "while the supplies last" constant-accretion scenario, 
we assume that the luminosity is a fraction £ of the Eddington 
luminosity LEddC^bhX an d that the luminosity is related to the 
mass accretion rate via L = eMc 2 , where as usual e parame- 
terizes radiative efficiency. The duration of the active phase 
is then Af act i ve = efm^c 2 / £L^i, implying a model luminosity 
iboi(M)h,0 = ^Edd#(Af ac tive - f X where H(x) is the Heaviside 
step function. 

2.3. Integration over the Black Hole Mass Function 

The LF of tidal disruption events is obtained by weighting 
the disruption probability modeled as a Poisson process by the 
black hole mass function. In the fallback model, the LF reads 

*(Lx)= / dM hh / dt *(M bh ) r(M bh ) e - r(M ^' 

x«y[w -1 Lboi(Afbh,0-^x], (4) 

where u> is the bolometric correction, 5(x) is the Dirac 6- 
function, and the dependence on / and to* is implicitly as- 
sumed; for the radius of sub solar stars we adopt the rela - 
tion R+ = Rq^/Mq) 8 (e.g.. lKippenhahn & Weigerl 199(H) . 
Here, M m \„ is the minimum MBH mass, and as before M max oc 
m°' 7 is the maximum MBH mass near which a star can be dis- 
rupted. We take M m m = 10 4 M Q , as systems in which an even 
smaller black hole is exp ected have been found to lack one 
(e.g., iMerritt et al.ll200lt iVallnri et aU fcoO^. For the bolo- 
metric correction in the 0.5 - 2 keV band, we take uj = 45, 
which at lo w lu minosities is approximately independent of 
luminosity ( Marconi et all 120041 and references therein), al- 
though in reality the bolometric correction could be correlated 
with Lboi/^Edd rather than L bo i. 

In Figure^ we plot the resulting LF for our fiducial model 
with a light curve that decays as a power-law L oc t~ x with 
A = |. The model is based on a galaxy LF with amplitude $0 = 
0.02 h~ 3 Mpc~ 3 and absol ute blue magn itude of an L+ galaxy 
of Mb* = -20.0 + 51og 10 /z jBlanton et afll2001l) . where we as- 
sume h = 0.7. We vary the various parameters, including A in 
Figure to explore the sensitivity of the X-ray AGN LF to 
our choice of the fiducial model. In the low-luminosity regime 
L x < Lpeak = ^ -I £boi[>w*,fpeak0ra*)], the LF is a power-law 
$(L X ) oc L x l/X ~ l . At L x > L peak ~ 10 44 (w/45)- 1 ml /3 ergs s" 1 , 
the LF drops precipitously. 

Below L x ~ 10 44 ergs s" 1 , the calculated X-ray AGN LF 
slightly overestimate s but is roughly com patible with the soft 
X-ray AGN LF of fHasinger et al. (2005) (shown in the fig- 
ure) ,andJis_afecjx>rof aTew^ hard X-ray AGN 
LF (lUeda et all 120031 iBarger et al]l2005l) . The discrepancy 
between the soft and hard X-ray LF could be attributed to an 



incompleteness of the soft X-ray LF resulting from obscura- 
tion. 

In models with power-law decaying light curves, the low- 
luminosity slope of the AGN LF is determined by the power- 
law of the decay, rather than the shape of the black hole mass 
function (Fig.Q/). The amplitude of the LF, however, is sen- 
sitive to the black hole mass function. The steep drop in the 
LF around 10 44 ergs s" 1 is associated with the peak luminos- 
ity produced in a tidal disruption event, which corresponds to 
the peak mass inflow rate from the debris that is beginning to 
circularize around the black hole (Fig.JQj,/?). 

In the "while the supplies last" model, where the accretion 
luminosity is a constant and equal to £Z-Edd, the LF reads 

*( Lx > ~ 1 T~7^~ 7^r(Mi)^ r (Mi), (5) 
1 4irum p Lx 

for M min < Mi < M max , where Mi = uHnL x /4ir£Gm p c. In 
the low luminosity regime, the LF is a power law ^(Lx) oc 
^+1 /2a+(a+\)k j/i-i /2a-(a+\)k ^ wn i cn assuming a constant £ is 
steeper, Lx^(Lx) oc L x 15 , than the observed relation (fi > 
— 1), although the two LFs have comparable amplitudes at 
Lx ~ 10 43 ergs s" 1 . The discrepancy may reflect a luminosity- 
dependent Eddington ratio £. The dependence of £ on lumi- 
nosity or black hole mass is expected if the accretion rate is 
supply limited, rather than radiatively limited. For example, 
to reproduce ^{Lx) oc L x 2 , we must have £ oc Lx- At present 
it is not clear how such dependence should arise. 

3. DISCUSSION 

While many theoretical uncertainties frustrate an accurate 
estimate of the contribution of activity associated with stellar 
tidal disruption to the X-ray LF of active galaxies, it seems 
that the contribution can be significant at low luminosities 
L x < 10 44 ergs s . The tidal disruption activity makes a neg- 
ligible contribution at higher luminosities (power-law profile 
of the measured soft X-ray LF extends to Lx ~ 10 47 5 ergs s" 1 
at redshifts z ~ 1-3). MBHs with masses above ~ 10 7 Mq 
grow by accreting non-stellar material from the interstellar 
medium. 

What is the origin of the featureless knee in the LF at 
Lx ~ (10 43 - 10 44 ) ergs s" 1 ? The location of the knee co- 
incides with Eddington-limited accretion onto black holes 
with masses M b h ~ 10 6 ~ 7 Mq.' Given the tidal disruption 
rate assumed here (eq. black holes with masses M b h < 
2 x 1O 6 (//O.25)M , where / is the fraction of the stellar mass 
that ultimately gets accreted onto the MBH, could have grown 
to their present size by accreting stellar debris over a Hubble 
time, although undoubtedly, non-stellar accretion must have 
played a role. This critical mass is compatible with the loca- 
tion of the knee in the LF assuming near-Eddington accretion. 

Black holes with masses ~ 1O 7 M , in turn, roughly lie at 
the demarcation line separating the dominance of late-type 
and early-type galaxies. If major mergers are responsible for 
diverting gas into accretion onto black holes in early type 
galaxies, then the absence of merging activity in late-type 
galaxies leaves stellar tidal disruption as competitive supplier 
of material to the MBH. This may explain why the steep 
power-law behavior of the quasar LF fails to extend to lower 
luminosities. 

The observed X-ray LF of AGN evolves with redshift. The 
number density of low-luminosity AGN increases by a factor 

1 The maximum accretion rate might exceed the Eddington limit by a fac- 
tor of a few (e.g., Beeelman 2002). 
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FIG. 1. — The X-ray LF of tidal disruptions calculated using equation |4} assuming a bolometric correction ui = 45 and a light curve that decays 
as a power-law A = |. In each panel, a parameter is varied to illustrate dependence: (a) m„ = (0. 1,0. 333,7.0, 3. 0)Mq; (b) a = (—1,-1.15,-7 .3, -1.45); 

(c) k = (0.70, 0.75, 0.8, 0.85); (d) A = 1.222,1.444,7.666,1.888; (e) M s * = (-18,-19,-20,-21) + 51og 10 /i; (f) M rain = (10 3 ,7O*, 10 5 , 10 6 )M©; (g) a = 
(4.0,4.333,4.666,5.0); and (h) £ = 0.5, 7, 2,4, where the fiducial value has been italicized. We also show the low-redshift (0.015 < z < 0.2) soft X-ray (0.5-2 keV) 
LF data from the Chandra and XMM-Newton survey of Hasinger et al. 1 2005). A Hubble constant of 77q = 70 km s -1 Mpc~' is assumed throughout. 



of - (10,4) at L x ~ (10 42 , 10 43 ) ergs s" 1 fromz ~0to z~ 0.7 
in the soft X-rays (0.5-2keV; Hasin ger et afll2005l) . while 
the respective incr ease is only a facto r ~ 2 - 3 in the hard 
X-rays (2-8 keV; iBarger et al J 120051) . The location of the 
knee shifts fromLx ~ 10 43 ergs s" 1 to Lx ~ 10 44 ergs s" 1 from 
Z ~ to z ~ 1 in the soft and hard band alike, and the break 
becomes more prominent at higher redshift. These trends re- 
quire explanation. 

The number density of tidal-disruption powered sources 
could increase with redshift because of one or more of the 
following reasons: 1. On average, MBHs are smaller at larger 
redshift, implying a larger tidal disruption rate or, specula- 
tively, a higher radiative efficiency in the soft X-ray band; 2 2. 
Galactic nuclei were denser at larger redshift (or, equivalently, 
stellar velocity dispersions were higher at fixed MBH mass), 
also implying a larger tidal disruption rate; and 3. Tidally 
disrupted stars were more massive at larger redshift, implying 
brighter events. 3 Any one or all of these factor may contribute 
to the observed trend at low luminosities, and could also be 
responsible for the apparent brightening of the location of the 
knee of the LF. 

In particular, if MBH grew gradually in pre-existing stellar 
density cusps with fixed velocity dispersions, then the tidal 
disruption rate would be inversely proportional to the black 
hole mass (eq.[0 and would increase with redshift. This could 
explain the increase in the number density of low-luminosity 
AGN from the nearby universe to z ~ 0.7. 

2 Note that for black holes smaller than < 10 6 Mq, thermal emission 
from an Eddington-limited disk enters the soft X-ray band, and the soft X-ray 
luminosity is thus anticorrelated with mass. It is also possible that transition 
from high to low-state accretion is sensitive to My, . 

3 The remaining possibility, an increase of the number of black holes in 



The relative contribution of stellar tidal disruptions to the 
X-ray AGN LF can be studied by examining long-term X- 
ray variability of the sources. Currently, long-term monitor- 
ing d ata have been analyzed for only a handful of sources 
(e.g., Markowitz & Edelson 2004). Theoretical uncertainty 
notwithstanding, the light curves of tidal disruption events 
should be characterized by a fast rise and slow decay, but 
the variability may or may not be detected as an outburst. 
The frequency of hard X-ray outburs t s in th e ROSAT All-Sky 
Survey was studied by Donl ev et al.l (|2£)02). Their estimate, 
~ 10~ 5 yr" 1 per galaxy, is significantly below the tidal disrup- 
tion rate adopted by us. The rates can be reconciled if the 
majority of tidal disruptions do not exhibit outburst character- 
istics. In that case, a measurable asymmetry in the autocorre- 
lation function will still be expected. 
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the centers of galaxies with redshift, seems implausible. 



Ayal, S., Livio, M., & Piran, T. 2000, ApJ, 545, 772 



REFERENCES 

Bahcall, J. N, & Wolf, R. A. 1976, ApJ, 209, 214 



MILOSAVLJEVIC ET AL. 



5 



Bahcall, J. N., & Wolf, R. A. 1977, ApJ, 216, 883 

Barger, A. J., Cowie, L. L., Mushotzky, R. F., Yang, Y., Wang, W.-H., Steffen, 

A. T„ & Capak, P. 2005, AJ, 129, 578 
Begelman, M. C. 2002, ApJ, 568, L97 
Blanton, M. R., et al. 2001, AJ, 121, 2358 
Cannizzo, J. K., Lee, H. M., & Goodman, J. 1990, ApJ, 351, 38 
Cohn, H„ & Kulsrud, R. M. 1978, ApJ, 226, 1087 

Donley, J. L„ Brandt, W. N., Eracleous, M., & Boiler, T. 2002, AJ, 124, 1308 

Elvis, M„ et al. 1994, ApJS, 95, 1 

Evans, C. R., & Kochanek, C. S. 1989, ApJ, 346, L13 

Faber, S. M., et al. 1997, AJ, 1 14, 1771 

Ferrarese, L. 2002, Current high-energy emission around black holes, 3 
Ferrarese, L., et al. 2006, in press 

Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Frank, J., & Rees, M. J. 1976, MNRAS, 176, 633 
Gebhardt, K., et al. 1996, AJ, 112, 105 
Genzel, R., et al. 2003, ApJ, 594, 812 

Halpern, J. P., Gezari, S., & Komossa, S. 2004, ApJ, 604, 572 
Hasinger, G, Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Hills, J. G. 1975, Nature, 254, 295 

Hopkins, P. E, Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, B., & 

Springel, V. 2005, ApJ, 630, 716 
Kim, S. S., Park, M.-G, & Lee, H. M. 1999, ApJ, 519, 647 
Kippenhahn, R., & Weigert, A. 1990, Stellar Structure and Evolution, (Berlin: 

Springer) 

Kobayashi, S., Laguna, P., Phinney, E. S., & Meszaros, P. 2004, ApJ, 615, 
855 

Kochanek, C. S. 1994, ApJ, 422, 508 
Khokhlov, A., & Melia, F. 1996, ApJ, 457, L61 
Komossa, S., & Greiner, J. 1999, A&A, 349, L45 
Komossa, S. 2004, IAU Symp., 222, 45 

Komossa, S., Halpern, J., Schartel, N., Hasinger, G, Santos-Lleo, M., & 
Predehl, P. 2004, ApJ, 603, L17 



Kroupa, P. 2001, MNRAS, 322, 231 
Lee, H. M., Kang, H., & Ryu, D. 1996, ApJ, 464, 131 
Li, L.-X., Narayan, R., & Menou, K. 2002, ApJ, 576, 753 
Lightman, A. P., & Shapiro, S. L. 1977, ApJ, 211, 244 
Magorrian, J., & Tremaine, S. 1999, MNRAS, 309, 447 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Marconi, A., Risaliti, G, Gilli, R., Hunt, L. K, Maiolino, R., & Salvati, M. 

2004, MNRAS, 351, 169 
Markowitz, A., & Edelson, R. 2004, ApJ, 617, 939 
Menou, K., & Quataert, E. 2001, ApJ, 562, L137 
Merritt, D., Ferrarese, L., & Joseph, C. L. 2001, Science, 293, 1 1 16 
Merritt, D., & Szell, A. 2005, preprint I astro-ph/05 10498 i 
Merritt, D., & Poon, M. Y. 2004, ApJ, 606, 788 
Meszaros, P., & Silk, J. 1977, A&A, 55, 289 
Milosavljevic, M., & Loeb, A. 2004, ApJ, 604, L45 

Nakamura, O., Fukugita, M., Yasuda, N., Loveday, J., Brinkmann, J., 

Schneider, D. P., Shimasaku, K., & SubbaRao, M. 2003, AJ, 125, 1682 
Phinney, E. S. 1989, IAU Symp. 136: The Center of the Galaxy, 136, 543 
Rees, M. J. 1988, Nature, 333, 523 
Rees, M. J. 1990, Science, 247, 817 

Schodel, R., Eckart, A., Genzel, R., Merritt, D., Alexander, T., Ott, T., 

Straubmeier, C. 2006, submitted 
Simien, E, & de Vaucouleurs, G. 1986, ApJ, 302, 564 
Ueda, Y, Akiyama, M., Ohta, K., & Miyaji, T. 2003, ApJ, 598, 886 
Ulmer, A. 1999, ApJ, 514, 180 

Valluri, M., Ferrarese, L., Merritt, D., & Joseph, C. L. 2005, ApJ, 628, 137 

Wang, J., & Merritt, D. 2004, ApJ, 600, 149 

Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 

Young, P. J., Shields, G. A., & Wheeler, J. C. 1977, ApJ, 212, 367 



